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Abstract. Lattice refinement in LQC, its meaning and its necessity are discussed. The role 
of lattice refinement for the realisation of a successful inflationary model is explicitly shown. A 
, simple and effective numerical technique to solve the constraint equation for any choice of lattice 

■ refinement model is briefly illustrated. Phenomenological and consistency requirements leading 

to a particular choice of lattice refinement model are presented, while it is subsequently proved 
that only this choice of lattice refinement leads to a unique factor ordering in the Wheeler-De 
Witt equation, which is the continuum limit of LQC. 



1. Introduction 

Loop Quantum Gravity (LQG) [TJ [2] is one of the most promising candidate theories for 
l/^ | describing the quantum degrees of freedom of the gravitational field. Quantum gravity, 

combining consistently quantum mechanics and general relativity, is essential when curvature 
becomes large, as for example in the early stages of the evolution of our universe. LQG is 
a non-perturbative and background independent canonical quantisation of General Relativity 
(GR) in four space-time dimensions. Considering a fundamental theory of quantum gravity 
and following the main concept of GR, namely that gravity is indeed geometry, it follows that 
there is no background metric, there is only a manifold, while geometry and matter should both 
have a quantum mechanical origin. This indeed differentiates LQG from other approaches of 
quantising gravity, which have been developed in the framework of particle physical While the 
full theory is not yet complete, LQG already has a number of successes, such as the construction 
of quantum geometry the prediction of Planck discreteness in geometric operators, as well as a 
quantum accounting for black hole entropy on the horizon. 

The application of ideas and mathematical methods of the full LQG theory to the cosmological 
sector — the dynamical variables are reduced to first homogeneity and then possibly also to 
isotropic models — led to Loop Quantum Cosmology (LQC) [3j, which is not a field theorjQ. 
LQC, which gains a constantly increasing interest from the scientific community and has recently 
made significant progress, differs from other attempts of quantising gravity in the sense that it 

1 Background independence means that quatisation is achieved in the absence of a metric other that the physical 
one determined by the densitised triad. At the classical level, the fact that the laws of physics are background 
independent is mathematically expressed by the Einstein equations being four-diffeomorphism covariant. 

2 The only known consistent perturbative approach to quantum gravity is string theory, a theory aiming at 
unifying all interactions. 

3 The quantised Hamiltonian is given in terms of symmetry reduced variables, so there is only a finite number of 
degrees of freedom. 



gets its input from the full theory of LQG. In short, it is a cosmological mini-superspace model 
quantised with methods of the full LQG theory. Considering a flat isotropic model within LQC, 
the extrinsic curvature scale k = a = y / 87rGa 2 A/3 (a stands for the scale factor and A denotes 
a positive cosmological constant) appears in the holonomies in such a way that only e iak , with 
a G R, can be represented as operators, k itself cannot The parameter a is related to the 
edge length used in holonomies, which are playing the role of the basic operators and they imply 
that the Hamiltonian constraint — a Hamiltonian density which is constrained to vanish by the 
equations of motion — is quantised to a second order difference equation, instead of the second 
order differential equation of the Wheeler-De Witt (WDW) approach to quantum cosmology. 

Following the approach of LQC is useful in two ways. Firstly, it allows us to get some useful 
insight about open issues of the full LQG theory, and secondly, by using symmetry reduction 
of the infinite dimensional phase space of LQG, the theory becomes often tractable leading to 
satisfactory answers about various interesting physical questions. The discreteness of spatial 
geometry, a key element of the full theory, leads to successes in LQC which do not hold on the 
WDW approach. In particular, it has been shown that classical big bang [5] and black hole 
singularities [6] are removed in LQC, in a well defined manner. 

In LQC, the evolution of the universe is divided into three distinct phases depending on the 
value of the scales probed by the universe itself. They are the following: 

• The discrete quantum phase: Very close to the Planck scale, the concept of space-time 
has no meaning and one should consider the full quantum gravity theory. Applying LQC 
during this phase, leads to a finite bounded spectrum for eigenvalues of inverse powers of 
the three-volume density, called the geometrical density. 

• The intermediate phase: As the volume of the universe increases with time, the universe 
enters a semi-classical phase. More precisely, for lengths above Lp\ = ■s/jlpi (7 ~ 0.2375 
is the Barbero-Immirzi parameter — a constant ambiguity parameter, whose value is fixed 
by black hole entropy calculations — and lp\ is the Planck length with = (8itG)~ 1 ) 
the space-time can be approximated by a continuous manifold and the equations of motion 
take a continuous form. This intermediate phase is characterised by a second scale L* with 
L+ = \J (7J/in)/3/pi, below which the geometrical density is significantly different from its 
classical forrdj- For scales below L* quantum corrections cannot be neglected. This phase 
is the most relevant one regarding phenomenological consequences of LQC. 

• The classical phase: At later times, and therefore large scales, the universe enters the full 
classical phase and standard cosmology becomes valid. 

LQG/LQC is formulated in terms of SU(2) holonomies of the connection and triads^ |7J. 
To obtain a quantum constraint we introduce an operator representing the curvature of the 
gravitational connection. A feature of LQC, which is a direct consequence of LQG, is that 
while there are well-defined analogs of holonomies, there is no operator corresponding to the 
connection; one defines a curvature operator in terms of holonomies. In the classical level, 
curvature can be expressed as a limit of the holonomies around a loop as the area enclosed by 
the loop shrinks to zero. However, in quantum geometry one cannot continuously shrink a loop to 
zero area, since the eigenvalues of the area operator are discrete, implying that there is a smallest 
non-zero area eigenvalue, called the area gap [3 [9]. In a canonical quantisation scheme, as the one 
followed in this approach, one first writes the action of GR into a Hamiltonian formulation and 

4 The half-integer J labels the ambiguity in choosing the representation in which the matter part of the 
Hamiltonian constraint for a scalar field is quantised. The length parameter fj,o is related to the underlying 
discrete structure. 

5 The connection in LQG determines the parallel transport of chiral fermions, mathematically represented by 
spinors, while the conjugate momenta can be interpreted as spatial triads (i.e., square-roots of the metric of the 
3-dimensional space). 



then one quantises this classical Hamiltonian. In the old quantisation, the quantised holonomies 
were taken to be shift operators with a fixed magnitude. However, it was shown |10|. [TT] that 
this approach leads to unavoidable instabilities in the continuum semi-classical limit, where the 
WDW wave-function becomes a good approximation to the difference equation of LQC. For a 
large semi-classical universe, the WDW wave-function would be oscillating on scales of the order 
(av^Vr 1 . As the universe expands, this scale becomes eventually smaller than the discreteness 
scale of the difference equation of LQC, implying that discreteness of spatial geometry would 
become apparent in the behaviour of wave-functions describing a classical universe. In the 
underlying LQG theory, the contributions to the discrete Hamiltonian operator depend on the 
state which describes the universe. As the universe expands, the full Hamiltonian constraint 
operator creates new vertices of a lattice state, in addition to changing their edge labels. As 
the extrinsic curvature scale k increases with increasing volume, the corresponding a decreases 
since the lattice is being continuously refined. Thus, in the context of LQC one has a refinement 
of the discrete lattic^f]. One can choose such a lattice refinement model, so that the increase 
in extrinsic curvature scale k can be balanced by the decrease of a such that ak remains small 
and semi-classical behaviour is achieved for any macroscopic volume even in the presence of a 
positive cosmological constant. 

Lattice refinement is also required from phenomenological reasons [111 I12j. For example, 
as we will discuss, lattice refinement renders a successful inflationary er4j more natural |12j . 
The effect of lattice refinement has been modelled and the elimination of the instabilities in the 
continuum era has been explicitly shown. 

The correct lattice refinement model should be obtained from the full LQG theory. In 
principle, one should use the full Hamiltonian constraint and find the way that its action balances 
the creation of new vertices as the volume increases. Instead, phenomenological arguments have 
been used, where the choice of the lattice refinement model is constrained by the form of the 
matter Hamiltonian [T7j. In particular, we have shown [T7] that for LQC to generically support 
inflation, and other matter fields, without the onset of large scale quantum gravity corrections, 
one should adopt a particular model of lattice refinement. This choice has been then found [18] 
to be the only one, for which physical quantities are independent of the choice of the elementary 
cell used to regulate the spatial integrations. Amazingly enough, this is exactly the choice 
required for the uniqueness of the factor ordering of the Wheeler-De Witt equation [19] . 

Lattice refinement leads to new dynamical difference equations which, in general, do not 
have a uniform step-size making their study quite involved, particularly for anisotropic cases, 
as for example for Bianchi models or black hole interiors. Numerical techniques have been 
developed [201 EI] to address this issue. 

Our primary interest here concerns cosmological predictions of quantum gravity, and therefore 
we only focus on questions dealing with the early universe and the initial conditions. A quantum 
theory of gravity is expected, on the one hand, to cure the classical singularities of GR and, on 
the other hand, either to provide the conditions suitable for the onset of inflation, or to suggest 
an alternative scenario for alleviating the strong fine-tuning of the standard cosmological model. 

6 The parameter, denoted earlier by no, which appears in the regularisation of the Hamiltonian constraint will 
no longer be constant within the lattice refinement context. 

7 It was hoped that LQC would help to overcome the extreme fine tuning necessary to achieve successful inflation 
in GR [13H14I IT5] . However, it has been shown [16] that semi-classical corrections are insufficient to alleviate this 
difficulty. Certainly, inflation could be generic in the deep quantum regime. 



2. Elements of LQG/LQC 

LQG, as well as LQC, are both based on a Hamiltonian formulation of GR with basic variables 
an SU(2) valued connection and the conjugate momentum variable which is a densitised triacd, a 
derivative operator quantised in the full LQG theory in the form of fluxes. By using connection- 
triad variables, arising from a canonical transformation of Arnowitt-Desner-Misner (ADM) 
variables, we make an analogy with gauge theories; this will be helpful when dealing with 
quantisation issues. 

The densitised triad carries information about the spatial geometry, encoded in the three- 
metric; the connection carries information about the spatial curvature, in the form of the spin- 
connection, and the extrinsic curvature. More precisely, the densitised triad, Ef, is related to 
the three- metric, q a b, by 

E? = \/W\e? , (1) 
where e? is a physical triad, dual {e\e° a = 6f) to the co-triad, ei, and satisfying 

lab = e* egdij . (2) 

Note that i refers to the Lie algebra index and a is a spatial index with a, i = 1, 2, 3. 
The connection, A l a , can be related to the ADM variables as 

Ai = Ti + 7 Ki, (3) 

where r a = -(l/2)e ijk e b {2d [a e k ] + e c k e l a d c e l b ) stands for the spin-connection compatible with the 
co-triad, 7 is the Barbero-Immirzi parameter — classically it has no physical consequence, while 
at the quantum level it plays a role in the level spacing of discrete geometric eigenvalues — and 
K l a stands for the extrinsic curvature one- form K l a = e bl K a b (with K a b the extrinsic curvature), it 
is the Lie derivative of q with respect to the normal vector to the spatial slice, K l a = (Cftq a b)S^ e b . 
The Poisson bracket of the densitised triad and the connection reads 

{A\{x) 1 E b {y)} = K 1 5 b a 5)5\x,y) , (4) 

where k = 8ttG. 

The Hamiltonian for GR is given by the sum of constraints, with the scalar constraint 

C GR = r 1 - / d 3 xiV(x) 
2k Jt, 

being the most important one. Note that i 7 ^ = 2<9[ a A^ + e^A^A-L denotes the curvature 
two-form of the connection, and N stands for the lapse function. 

Before proceeding, let us note that for any quantisation scheme based on a Hamiltonian 
framework, as for example LQC, or an action principle, as for example path integral approaches, 
for the homogeneous flat {k = 0) model, one should regularise the divergences which appear due 
to the homogeneity as the action and Hamiltonian are integrated over spatial hyper-surfaces. 
To do so, the spatial homogeneity and Hamiltonian are restricted to a fiducial celO, with finite 
volume Vq = J d 3 Xy/pg[ [1]. 

8 The variables i?" and A x a were introduced by Barbero [22] as an alternative to the complex Ashtekar [23] 
variables. Both real and complex connections have been used for canonical gravity. The complex connection has 
SL(2,C) as gauge group, while the real connection has SU(2) as gauge group. Mathematical techniques can only 
cope with a quantum theory based on SU(2). 

9 Note that /10 , which we will discuss later is the scale of the finite fiducial cell that spatial integration is restricted 
to, so as to remove the divergences that occur in non-compact topologies. 
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The expressions can be simplified a lot by restricting the analysis to homogeneous and 
isotropic geometries: 

q = [a{t)] 2 5 l3 °oji °4dx a dx b , (6) 
where the one-forms satisfy dy a = 0, leading to 

q = [a(t)] 2 [(l - kr 2 )- 1 ^ 2 + r 2 d0 2 ] . (7) 

The symmetry reduced variables E, A are given by 

A\ = V~ 1/3 c °u4 , E? = pV~ 2/3 \det °u>\ °e? , (8) 

where the vector fields °e? are dual to the 1-forms: °u)l °e? = 5\. 

Concentrating on isotropic cosmological models with constant spatial curvature, we will first 
consider symmetric connections and triads, and we will then insert them into the full action 
leading to a symmetry reduced action. The symmetric connections and triads can be decomposed 
using basis one- forms and vector fields obtained by Bianchi models. One can show that the 
actions lead to the correct equations of motion of GR, implying that they are equivalent to the 
Einstein-Hilbert action on metric variables. 

The loop quantisation of the flat FLRW scalar constraint changes the curvature 2-form F, 
its ab component reads [24] : 

F k ab = -2 lim^oTr ( -9^- | r k V 6 ; (9) 

Ar^\ is the area of the square dy in the (i,j)-plane swept by a face of the elementary cell, the 

holonomy /if? ' around the square Djj is the product of holonomies along the four edges of D^- , 

and p, is p = y/A/lpl 1 ^ 2 , with A the eigenvalue of the area operator. 

In the case of a spatially flat background, derived from the Bianchi I model, the isotropic 
connection can be expressed in terms of the dynamical component of the connection c{t) as 

4 = m4 > (io) 

with uj l a a basis of left-invariant one- forms uj l a = dx l . The densitised triad can be decomposed 
using the Bianchi I basis vector fields Xf = 6f as 



Ef = V°qp(t)x? , (n) 

where °q stands for the determinant of the fiducial background metric, 

°q ab = J a u u , (12) 

and p{t) denotes the remaining dynamical quantity after symmetry reduction. 

In terms of the metric variables with three-metric q a b = a 2 uj l a u>bi, the dynamical quantity is 
just the scale factor a(t). Given that the Bianchi I basis vectors are Xf = Sf, 

\p\ = a 2 , (13) 

where the absolute value is taken because the triad has an orientation. Since the basis vector 
fields are spatially constant in the spatially flat model, the connection component is 



c = sgn(p)7— 



(14) 



Note that in what follows, the lapse function which is a constant due to spatial homogeneity, 
will be set equal to 1. Thus, GR can be formulated as a gauge theory in Ashtekar variables. 
The canonical variables c,p are related through 

{c,P} = ^Vo , (15) 

where Vo the volume of the elementary cell adapted to the fiducial triad. 

Defining the triad component p, determining the physical volume of the fiducial cell, and 
the connection component c, determining the rate of change of the physical edge length of the 
fiducial cell, as 

p = V* /3 p , c = F 1/3 c, (16) 

respectively, we obtain 

{c, P } = ^, (17) 

independent of the volume Vo of the fiducial cell. 

To quantise, we follow the approach used in the full LQG theory. The metric itself is a 
physical field which must be quantised; it cannot be considered as a fixed background. Thus, 
to quantise gravity we use gauge theory variables to define holonomies of the connection along 
a given edge 

h e (A) =Vexpj d S 7 M 004W s )) r * > ( 18 ) 

where V indicates a path ordering of the exponential, 7^ is a vector tangent to the edge and 
Tj = —iai/2, with <7j the Pauli spin matrices, and fluxes of a triad along an S surface 

E(S, f)= [ e abc E ci f idx a dx b , (19) 

with fi an SU(2) valued test function. Note that even though these variables appear rather 
artificial, as presented here in the context of LQC, they nevertheless arise naturally within the 
full LQG theory. 

Thus, the basic configuration variables in LQC are holonomies of the connection 

^ 0) (^)=cos(^)l + 2sin(^)r,, (20) 

along a line segment /jlq °ef and the flux of the triad 

F s (E,f)<xp; 

the basic momentum variable is the triad component p. Note that 1 is the identity 2x2 matrix 
and Tj = —iai/2 is a basis in the Lie algebra SU(2) satisfying the relation 

T iTj = (l/2)e ijk r k - (l/4)<% . 

Let us first review the old quantisation procedure. We follow the same approach as in LQG. 
We thus take (with an arbitrary real number) and p, as the elementary classical 

variables, which have well-defined analogues [3J. Using the Dirac bra-ket notation and setting 
gi^oc/ 2 = (c|/z), the action of the operator p acting on the basis states \fi) is 



(21) 



where fi (a real number) stands for the eigenstates of p, satisfying the orthonormality relation 



<A*l|A*2> 



(22) 
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operator acting on basis states \fi) is 
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i Ho 



\H) = exp 



= |m + Mo) , 



(23) 



where [io is an Y rea l number. Thus, in the old quantisation, the operator e wc / 2 acts as a simple 
shift operator. 

Using the volume operator V = |j5[f^ 2 , representing the volume of the elementary cell with 
eigenvalues = (k7^|/j|/6) 3 / 2 , we ge£3 



3/2 



(24) 



To define the inverse volume operator one has to trace over SU(2) valued holonomies. Since 
there is a freedom in choosing the irreducible representation to perform the trace, an ambiguity 
J arisea^l Let us use the J = 1/2 irreducible representation of SU(2). The inverse volume 
operator is diagonal in the basis and is given by [I] 



6 



{V(» + »o)} 1/3 -{V(»-fio)} 1/3 



I A*) 



(25) 



where no is proportional to the length of the holonomy. Note that the regulating length [iQ is 
the crucial parameter in the quantum corrections; (jlq determines the step-size of the resulting 
difference equation. In the above equation, Eq. (|25p. the eigenvalues are bounded and approach 
zero near the classical singularity; in the classical case the eigenvalues diverge at the singularity 
[i = 0. The eigenvalues reach their maximum at a characteristic scale /j,q, at larger scales they 
approach the classical values and at smaller scales they are suppressed [25] . 

As in the full LQG theory, there is no operator corresponding to the connection. Nevertheless, 
the action of its holonomy is well-defined. Let us denote by tif" ^ the holonomy along the edge 



parallel to the i basis vector of length HoVq 
on the basis states is given by [26] 



1/3 



with respect to the fiducial metric. Its action 



^°V) = (&l + 2sn7i)|/i) 



(26) 



where, 



cs|/x) = cos(/x c/2)|/i) 
sh\n) = sin(/i c/2)|/i) 



[ l^ + ^o) + Ia* - 
[ \n + Ho) - Ia*- 



A*o; 
Mo) 



/2, 
/(2i) 



(27) 



Thus, 



^o)^o)^Mo)-l^o)-l|^ 

= [(cs 4 - sn 4 ) 1 + 2(1- ArjTi) cs 2 sn 2 + 4 (r^ - r») lcs 



sn 



(28) 



10 Being concerned with the large scale behaviour of the LQC equations, we neglect the sign ambiguity. 

11 Note that J/2 stands for the spin of the representation. Usually one quantises the gravitational part of the 
Hamiltonian constraint using the fundamental J — 1/2 representation, and the ambiguity is only investigated for 
the matter part. 



and 



;(w>) 
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V — cst^cs — snFsn ) + 2T* ( cst^sn — snFcs ) | /i) . 



(29) 



The gravitational part of the Hamiltonian operator in terms of SU(2) holonomies and the triad 
component, in the irreducible J = 1/2 representational , reads 

2i 



C. 



grav 



tr ^ (^h^hi^-'h^-'hi^ 



ijk 



h 



(mo)-i 



V 



(30) 



As in the full LQG theory, curvature is defined in LQC in terms of holonomies around closed 
loops. This implies that the limit /xo — > of the above operator does not exist, since it would 
mean that the area enclosed by loops should be shrunk to zero. In the underlying quantum 
geometry, the eigenvalues of the area operator are discrete, implying that there is a smallest 
nonzero eigenvalue, an area gap A [9]. This is indeed the reason for which the WDW differential 
equation gets replaced by a difference equation whose step size is controlled by A. Since no 
enters through the holonomies, its value in the fixed lattice case was fixed by demanding that 
the eigenvalue of the area operator be the area gap: A = 2v / 37T7^p 1 , implying /io = 3\/3/2. 

The action of the self-adjoint Hamiltonian constraint operator, 7Y gr av = (Cgrav + Cgrav)/2, on 
the basis states, is 

3 



Wgravl^) 



where 



4k 2 7 3 /i/Xq 



[R(jm) + R{n + 4/xo)] \n + 4mo) - 4R(ti) I ») + [R(j*) +R(n-4no)]\fi- 



4/Uo)} , 
(31) 
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(32) 



Having the Hamiltonian operator, dynamics are determined by the Hamiltonian constrain 
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grav 







(33) 



Note that in the full LQG theory, there is an infinite number of constraints, whereas in the 
reduced homogeneous and isotropic case there is only one integrated Hamiltonian constraint. 
Matter is then introduced by just adding the actions of matter components to the gravitational 
action. One finally obtains difference equations analogous to the differential WDW equations F^l 
Let us impose the constraint equation on the physical wave- functions \*&), which can be 
expanded using the basis states as |^) = ^V(^) !/•*)> with summation over values of fi and 
where the dependence of the coefficients on q> represents the matter degrees of freedom. Since 
the states \fj) are eigenstates of the triad operator, the coefficients $ ' n{4>) represent the state 
in the triad representation. Thus, quantising the Friedmann equation along the lines of the 
constraint in the full LQG theory, one gets the following difference equation |27j 
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HM*^) , (34) 



12 With this choice, the Hamiltonian constraint is free of the ill-behaving spurious solutions. 

13 Note that the Gauss and the diffeomorphism constraints are automatically satisfied by appropriate gauge fixing. 

14 The reader should note that in LQC the physical fundamental object is the discrete difference equation, while 
the differential equation is just the approximation in the continuum limit. 



where the matter Hamiltonian Ti.^ is assumed to act diagonally on the basis states with eigenvalue 
Ti^fi). Equation (fMJ) is indeed the quantum evolution (in internal time fi) equation. There 
is no continuous variable (the scale factor in classical cosmology), but a label [i with discrete 
steps. The wave- function ^^(<j)), depending on internal time ji and matter fields cj>, determines 
the dependence of matter fields on the evolution of the universe. A massless scalar field plays the 
role of the emergent time. Thus, in LQC the quantum evolution is governed by a second order 
difference equation, rather than the second order differential equation of the WDW quantum 
cosmology. As the universe becomes large and enters the semi-classical regime, the WDW 
differential equation becomes a very good approximation to the difference equation of LQC. 

3. Lattice refinement 

Consider the continuum limit (namely that fx S> liq) of the Hamiltonian constraint operator 
acting on the physical states. In the small regulating length fiQ limit, one obtains [17J a second 
order difference equation which distinguishes the components of the wave-functions in different 
lattices of spacing Alio. Assuming that VP does not vary much on scales of the order of A/xq, 
known as pre-classicality [28J, one can smoothly interpolate between the points on the discrete 
function ^(0) and approximate them by the continuous function ^(/x, eft). In this way, one 
approximates the difference equation by a differential equation for a continuous wave-function. 

The form of the wave-functions indicates that the period of oscillations can decrease as the 
scale increases, which implies that at sufficiently large scales the assumption of pre-classicality 
breaks down 15 !. This would then lead to quantum gravity corrections at large scales (i.e., 
classical) physics. To avoid this undesired event, was one of the motivations behind lattice 
refinement \27\ 129] . Allowing the length scale of the holonomies to vary, the form of the 
difference equation changes. Assuming that the lattice size is growing, the step-size of the 
difference equation is not constant in the original triad variables. The exact form of the difference 
equation depends on the lattice refinement used. 

Consider the particular model 



which is no longer a simple shift operator since fx is a function of fx. Changing the basis to 



Mo -> m(m) = MoM 



1/2 . 



(35) 



we will come back to the issue of the lattice refinement choice in a subsequent section. 
The basic operators are given by replacing liq with fx. Upon quantisation [2"fi] 



(36) 




(37) 



one gets 




,■- d 



d 



(38) 



The volume operator acts on these basis states as 




(39) 



The constant lattice, in the old quantisation approach, does not take into account the expansion of the fiducial 
cell in a FLRW background. 



and the self-adjoint Hamiltonian constraint operator acts as |25] 
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where 



C7 (i/) = |z> + /i | - \v - Mo I • 
Expanding |\&) = X^ 1 ^^)!^) the Hamiltonian constraint reads [T7] 
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Tiff, stands for the matter part of the Hamiltonian, which for a massive scalar field is given by 



(43) 



with Pff, the momentum and V((f>) the potential of the scalar field (p. Quantising the Hamiltonian 
constraint we obtain, in terms of p = 



dp 2 



dp2 (VP*(P, 0)) " 3p- 3/2 ^|^ + -^P^Vmtp, </>) + O (/i ) + • • • = 



_6_ 



which is a particular factor ordering of the WDW equation for a massive scalar field. 



(44) 



3. 1 . Implications for inflation 

We will show that lattice refinement is essential in order to achieve a successful inflationary era. 
Let us consider a fixed and a dynamically varying lattice and solve in the continuum limit the 
second order difference equation which governs the quantum evolution. By doing so we constrain 
the potential of a scalar field (the inflaton) so that the continuum approximation is valid. A 
second constraint is imposed on the inflationary potential so that there is consistency with 
measurements of the cosmic microwave background temperature anisotropics on large angular 
scales. Combining the two constraints in a particular inflationary model, for either a fixed lattice 
or a given lattice refinement model, we deduce the conditions for natural and successful inflation 
within LQC in each of the two cases. 

More precisely, let us separate the wave-function ^f(p,cj)) into Vl/(p, <f>) = Y(p)$(</>) and 
approximate the dynamics of the inflaton field, (f), by setting V{4>) = V^,p s ~ 3 ^ 2 , where Va, is 
a constant and 5 = 3/2 in the case of slow-roll, to get [17j 



V 



-1/2 _d_ 
dp 



p- 1/2 ^(p 3/2 T(p) 



+ pv^rip) = o 



(45) 



with solutions [T7] 



T (jp) p 
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7T 



I 37T 7T\ . / 37T 7T 
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(46) 
where 

x = 4v^(25 + 3 )-y 2<5 + 3 )/ 4 , (47) 

and = 96/ \nh 2 ). 

Without lattice refinement, the discrete nature of the underlying lattice would eventually 
be unable to support the oscillations and the assumption of pre-classicality will break down, 
implying that the discrete nature of space-time becomes significant on very large scales. For the 
end of inflation to be describable using classical GR, it must end before a scale, at which the 
assumption of pre-classicality breaks down and the semi-classical description is no longer valid, 
is reached. We will quantify this constraint. 

The separation between two successive zeros of T (p) is 

A P = -^=P {1 ~ 2S)/A ■ (48) 

For the continuum limit to be valid, the wave-function must vary slowly on scales of the order 
of n c = 4/2. Thus, we impose the constraint [T7] 



6~J 



Ap>4^^\ /+'-. (19) 
which implies the following constraint on Va, [Vi 



For slow-roll inflation, V(4>) must be approximately constant, implying 5 ~ 3/2. For /j,q = 3\/3/2 
and 7 rj 0.24, the constraint on the inflationary potential in units of h = 1 reads 

V(<j>) & 2.35 x Kr 2 Z pi 4 . (51) 

This is a softer constraint than the one imposed for fixed lattices, namely |17| 

« lO- 28 /-/ , (52) 

assuming that half of inflation takes place during the classical era. 

Selecting a successful and simple inflationary model, for example V(cf>) = m 2 </> 2 /2, an 
additional constraint can be imposed on the potential so that the fractional over-density in 
Fourier space at horizon crossing is consistent with the COBE-DMR measurements. Combining 
the two constraints we obtain [17] for the fixed and varying lattices 



m <, 70(e- 2iVcl ) Mpi (53) 
and m <, 10 M P1 , (54) 

respectively. 

Thus, for any significant proportion of inflation to take place during the classical era, the 
constraint imposed on the inflaton mass is very strong, while it becomes natural once lattice 
refinement is taken into account. 



3.2. Relation between lattice refinement model and matter Hamiltonian 

A particular lattice refinement model can support only certain types of matter [12]. To prove 
it, let us parametrise the lattice refinement by a parameter A and the matter Hamiltonian by 
a parameter 5 and solve the Hamiltonian constraint. The restrictions on the two-dimensional 
parameter space will become apparent once we impose some physical restrictions to the solutions 
of the wave- functions |12j . 
More precisely, assuming 

fi = VoV A , (55) 

we obtain 

~^~ A (56) 



Ml - A) 
The WDW constraint equation reads 



* = 0. (57) 



Being interested in the large scale limit, we can approximate [12\ the matter Hamiltonian with 
TifP = v s e ((f)), leading to 

e (0) = e (0) = - is- 5 H g r^ ■ (58) 

A necessary condition for the wave- functions to be physical is that the finite norm of the physical 
wave-functions, which is defined by J^ =( a dv\v\ s ^i 1 'S>2, must be independent of the choice of 
4> = (/>o • The solutions of the constraint are renormalisable provided they decay, on large scales, 
faster than v~W 5 \ 

To solve the constraint equation, we need to specify the from of TC^, which has in general two 
terms with different scale dependence. Being interested in the large scale limit, one will be the 
dominant term, allowing to write [12] 

f3H^ = e v (<j))v 5 » , (59) 

where the function e u is constant with respect to v. 

We then solve [12] the constraint equation and consider only those solutions which are the 
physical ones. The large scale behaviour of the wave-functions must be normalisable — a 
necessary condition for having physical wave-functions — and the wave- functions should preserve 
pre-classicality at large scales — a necessary condition for the validity of the continuum limit [12] . 
This procedure leads to constraints to the two-dimensional parameter space (A, 5), which we 
illustrate below [12] in Fig. 1 for A in the range < A < —1/2, imposed from full LQG theory 
considerations [TT] . 

We thus conclude that the continuum limit of the Hamiltonian constraint equation is sensitive 
to the choice of model and only a limited range of matter components can be supported within 
a particular choice |12j . 



3.3. Numerical methods for solving the constraint equation for any lattice refinement model 
The lattice refinement implies new dynamical difference equations, which are not expected to 
have a uniform step-size, leading to technical complications. This becomes apparent in the case 
of two-dimensional wave-functions, such as those necessary to study Bianchi models of black 
hole interiors. More precisely, the information needed to calculate the wave- function at a given 
lattice point is not provided by previous iterations. We prescribe below a method [21] based 
on Taylor expansion that can be used to perform this desired interpolations with a well-defined 
and predictable accuracj0. 

16 A simple local interpolation scheme to approximate the necessary data points, allowing direct numerical 
evolution of two-dimensional systems has been also proposed in Ref. [20] . 




Figure 1. The allowed types of matter content are significantly restricted. Note that in the 
case of a varying lattice (A ^ 0) it is not always possible to treat the large scale behaviour of 
the wave-functions perturbatively (dashed line with crosses) |12j . 



Let us first note that the Hamiltonian constraint for a one-dimensional difference equation 
defined on a varying lattice, can be mapped onto a fixed lattice simply by a change of basis [21 pi. 
However this method is of no help for the two-dimensional case, where the Hamiltonian constraint 
is a difference equation on a varying lattice [H] , 



+C {H, T) [{fl + 26J 1V+4V,t - 2 (1 + 2 7 2 J2) ^ T + (fj,— 2<g ^-4^,r] 

S 6^ 

+ C_ (/X, T) [*^26»,t-26 t ~ *»+25»,t-26 t ] = yH^ w , (60) 

with 

C± = 2<J /1 ( v /|r±2<y T | + v /M) , (61) 
C = Vk + «5r| - \/k-<*T| , (62) 

where we have defined <5„ and 5 r as the step-sizes along the [i and r directions, respectively. 
The parameter 5, with < 5 < 1, gives the fraction of a lattice edge that the underlying graph 
changing Hamiltonian uses [11] . 

In the case of lattice refining, 5^ and 5 T are decreasing functions of \i and r, respectively, 
and the data needed to calculate the value of the wave-function at a particular lattice site are 
not given by previous iterations. One can use Taylor expansions to calculate the necessary data 

17 We have checked |21j the validity of our Taylor expansion numerical method in calculating the wave-function 
by comparing our results with those obtained by mapping the one- dimensional difference equation defined on a 
varying lattice, onto a fixed lattice by performing a change of basis. 



points [21]. More precisely, let us assume that the matter Hamiltonian acts diagonally on the 
basis states of the wave-function, namely 



H. 



(63) 



/'■ T 



Given a function evaluated at three (noncolinear) coordinates, the Taylor approximation to the 
value at a fourth position is 
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(64) 



where the Taylor expansion is taken about the position (x2, 2/2)5 w ^ have defined 5f- 



x 4 



and Sfj = yi 



yj, and the differentials can be approximated using 



df 

f(xi,yi) = / (2:2,2/2) +6 i2q^ 
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(65) 
(66) 



where the dots indicate higher order terms. 

For slowly varying wave-functions, it has been shown [21] that linear approximation is very 



accurate and higher order terms in Taylor expansion can only improve the accuracy by 10 %. 
This method can be applied in any lattice refinement model, while its accuracy can be estimated. 
Even though we have illustrated it in the case of black hole interiors, this method can be applied 
to anisotropic Bianchi models and in general to systems with anisotropic symmetries. 

By using this Taylor expansions method, we were able to confirm [21] numerically the stability 
criterion of the Schwarzchild interior, which was earlier found [TTJ using a von Neumann analysis, 
and investigate [H] how lattice refinement can change the stability properties of the system. 
Finally, the underlying discreteness of space-time leads to a twist [21J in the wave-functions, for 
both a constant lattice, as well as lattice refinement models. 



3.4- Uniqueness of WD W factor ordering and the lattice refinement choice 
The correct lattice refinement model should in principle be given by the full LQG theory. In 
this sense, one should consider the full Hamiltonian constraint and find the way that its action 
balances the creation of new vertices while the volume increases. Instead of doing so, as we have 
already discussed earlier, phenomenological arguments [12], [T7] have been used to constrain the 
choice of the lattice refinement model by the form of the matter Hamiltonian. Later on, it has 
been argued [18] that only the lattice refinement model fi = \s,§yT x l 2 [19j . can be achieved by 
physical considerations of large scale physics and consistency of the quantisation structure. We 
will show below that this choice is also the only one which makes the factor order ambiguities 
of LQC to disappear in the continuum limit [19] . 

Indeed, there are many ways of writing the gravitational part of the Hamiltonian constraint in 
terms of the triad and the holonomies of the connection, our quantisable variables. Writing [26] 
for example, 

) , (67) 

'" "' " i)k 

one immediately realises that there are many possible choices of factor ordering that could have 
been made at this point, since classically the actions of the holonomies commute. However, each 



■-grav 



1 



of these factor ordering choices leads to a different factor ordering of the WDW equation in the 
continuum limit. 

The action of the factor ordering chosen in Eq. (|67p above gives |19| 



e ijk tr (hihjhi l h- ] 1 h k h^ 1 ,V ) 



-24sh 2 cs 2 ^csFsn — shycs 



(68) 



while other choices have different action, as it has been explicitly found in Ref. [19j . 

Defining V\v) = V v \v), the action of the above factor ordering on a general state in the 
Hilbert space given by |^) = Yli/^W) rea ds [19] 



e ijk tr f hihjhr 1 ^ 1 ^ h~ L ,V J 



-3i 
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+ (Vu+5k — V u+ 3kj 1py+ik 



(69) 



Similarly, one can find [19] the action of any other factor ordering choice. 

Then one can take the continuum limit of these expressions by expanding ip v ~ ip (y) as a 
Taylor expansion in small k/v. For the factor ordering chosen for this illustration here, the large 
scale continuum limit of the Hamiltonian constraint reads [191: 



lim ejjktr [h-Ji^ 1 h- x /ik 

k/v— >0 ' 
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(70) 



Setting a = 3/xo/(2fc) and hq = k, all lattice refinement models will lead to the same continuum 
limit for the WDW equation, only for A = —1/2 [T3], in which case the WDW equation reads [TU] 



lim C grav |^) 

k/u^0 
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d 2 4>. 
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S dv 2 
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(71) 



Thus, there is only one lattice refinement model, namely fl = ^ofi" 1 / 2 , with a non ambiguous 
continuum limhQ. 

In conclusion, phenomenological and consistency requirements lead to a particular lattice 
refinement model, implying that LQC predicts a unique factor ordering of the WDW equation 
in its continuum limit. Alternatively, demanding that factor ordering ambiguities disappear in 
LQC at the level of WDW equation leads to a unique choice for the lattice refinement model. 



4. Conclusions 

LQG canonically quantises space-time via triad and holonomies of the connection. Full 
understanding of the theory has not yet been reached, nevertheless symmetry reduction versions 
akin to WDW mini-superspace model have been successfully developed. 

18 The A — —1/2 choice can be easily understood with the following simple argument [30]. In LQC, the basis 
sates are \fi) and the physical area of the fiducial cell is yV^ 3 . Consider, in the full LQG theory, N fluxes passing 
through a side of the fiducial cell and divide its surface in A elementary surfaces. Then 7VA = [iV^ 3 , implying 
TV = fiVo /3 /A. The holonomies of the connection are e iAc/2 = e iA ^ 1/3 / 2 an d the fiducial area is A 2 = Vq /3 /N. 
Thus, A <x m~ 1/2 - 



As a first approximation the quantised holonomies were taken to be shift operators with 
a fixed magnitude. This results in the quantised Hamiltonian constraint being a difference 
equation with a constant interval between points on the lattice. These models lead to serious 
instabilities in the continuum semi-classical limit. 

In the underlying LQG theory, the contributions to the discrete Hamiltonian operator depend 
on the state which describes the universe. As the universe expands, the number of contributions 
increases, so that the Hamiltonian constraint operator is creating new vertices of a lattice state, 
leading to a refinement of the discrete lattice in LQC. 

The lattice refinement can be modelled and the instabilities in the continuum era get 
eliminated. We have discussed here why lattice refinement seems to be necessary to achieve a 
natural inflationary era, and we have illustrated that only a limited range of matter components 
can be supported within a particular lattice refinement choice. We have then shown that 
factor ordering ambiguities in the continuum limit of the gravitational part of the Hamiltonian 
constraint disappear only for a particular choice of lattice refinement. 

Whilst the continuum limit of the lattice refinement models can be taken, there is a 
complication in directly evolving two-dimensional wave-functions, such as those necessary to 
study Bianchi models or black hole interiors. The information needed to calculate the wave- 
function at a given lattice point is not provided by previous iterations. We have shown that 
Taylor expansions can be used to perform this interpolation with a well-defined and predictable 
accuracy. We have then discussed how lattice refinement can alter stability conditions of the 
system. 

We have only focused here in a few aspects of LQC concerning lattice refinement. There is 
certainly a much vaster arena of themes within LQC, an area which is gaining a lot of interest 
from the scientific community. Its basic advantage is that it allows us to successfully address 
some physical features of our universe, while it gives us some valuable insight for the full LQG 
theory. 
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